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Abstract 

We connect the quasinormal modes corresponding to Dirac fermions in various 
black holes backgrounds to an = 2 supersymmetric quantum mechanics algebra, 
which can be constructed from the radial part of the fermionic solutions of the Dirac 
equation. In the massless fermion case, the quasinormal modes are in bijective corre- 
spondence with the zero modes of the fermionic system and this results to unbroken 
supersymmetry. The massive case is more evolved but as we prove, supersymmetry 
remains unbroken even in this case. 

Introduction 

Black holes in equilibrium are generally speaking, simple objects due to the very few pa- 
rameters that are needed to describe them [IH^- However it is physically impossible to 
have an isolated black hole in nature, due to the fact that matter always exists around 
them, interacting directly with the black hole. This leads us to the conclusion that a black 
hole is always in a some sort of perturbed state, with more parameters to describe it, 
than in the unperturbed state. The perturbation of a black hole results to the emission of 
gravitational waves. Quasinormal modes [TH32] describe a long period of damped proper 
oscillations of gravitational waves, and are extremely important in various physical phe- 
nomena. The direct observation of black holes is actually based on quasinormal modes, 
with the most dominating one being the fundamental mode, namely the lowest frequency 
in the spectrum. Supermassive, intermediate mass and stellar mass black holes are essen- 
tial to the observation of astrophysical phenomena. In four dimensions, the only solution 
of the Einstein equations, that describes a stationary rotational black hole is the Kerr met- 
ric. The quasinormal modes are the characteristic sound of black holes and being such, 
matter field perturbations of such extreme gravitational backgrounds can reveal many 
important properties of black holes. In addition they depend only on the mass, angular 
momentum and charge of the black holes thus enabling us to directly identify the spacetime 
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parameters. The field perturbation of a black hole consists of three time evolution stages, 
namely the wave burst originating directly from the source, the quasinormal mode stage 
and the power-law tail [27]. Research in the area of quasinormal modes of black-holes has 
become quite scientifically compelling the last decade [28]. This is due to the vast number 
of applications and implications of many theoretical frameworks that incorporate quasi- 
normal modes. Firstly the existence of a black hole can be directly verified by observing 
it's fundamental quasinormal mode. Additionally the thermodynamic properties of loop 
quantum gravity (an appealing alternative to string theory) black holes can be further 
understood using the quasinormal modes (the imaginary part of the quasinormal modes 
of the scalar field perturbation is equal to the Barbero-Immirzi parameter [26], see and 
references therein) . Moreover the quasinormal modes of anti de Sitter black holes have a 
dual physical correspondence to quantities of the dual conformal field theory via the well 
known AdS/CFT correspondence [33j. From astrophysical aspects, the most interesting 
spacetimes are the asymptotically flat ones, however the observation of the universe's ex- 
pansion motivated the study of quasinormal modes in de Sitter and anti de Sitter space 
times [30]. Quasinormal modes can yield which gravitational systems are stable under 
dynamical perturbations. Actually a static or non-static solution describing a compact 
object is stable if all it's quasinormal modes are decaying in time, on the contrary even if 
one mode is growing, the gravitational system is unstable [1]. 

The perturbation of a black hole can be achieved either by directly perturbing the grav- 
itational background or by simply adding matter or gauge fields in the black hole space- 
time [1] . In this paper we shall use the later approach, in the linear approximation, which 
suggests that the field has no back-reaction on the metric. Particularly, we shall study 
Dirac fermion systems around various black hole environments and study when the sys- 
tem possesses an unbroken hidden supersymmetry. Supersymmetry has been connected 
to quasinormal modes spectra in the past, but in a different context |34H38| . Most of 
these works studied bosonic quasinormal modes and their relation to supersymmetry. We 
study the zero modes of the fermionic system and directly relate these to the quasinormal 
modes. As we shall see, the zero modes and quasinormal modes have a bijective correspon- 
dence, a fact that can actually be very crucial for supersymmetry to be unbroken. The 
specific type of supersymmetry that we found to be related to the fermion system is an 
N = 2 super symmetric quantum mechanics with zero central charge. The aforementioned 
supersymmetry is inherent to many gravitational systems [39^40] . The problem of calcu- 
lating the fermionic quasinormal modes can be very difficult. Many well known numerical 
techniques can be used [il [25ti32| . It is therefore useful to have a complete picture of 
the quasinormal spectrum and related issues regarding Dirac spectral problems in various 
gravitational backgrounds. This paper studies supersymmetries of extreme gravitational 
backgrounds and relates the most possible observable quantities, the quasinormal modes, 
to the corresponding Witten index of the supersymmetric algebra. 

This paper is organized as follows. In section 1, we present in brief the most necessary 
tools of the N = 2 supersymmetry algebra and related issues, that will be needed in the 
rest of the paper. In section 2, we study a massless Dirac fermion in a Kerr gravitational 
background. Since the Kerr solution is the most physically interesting black hole satisfying 
the Einstein's equations we shall study this in great detail. In section 3, we study the mass- 
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less and massive Dirac fermion in Kerr-Newman, Reisnerr-Nordstrom and Schwarzschild 
gravitational backgrounds. The massive case proves to be much more evolved compared 
to the massless case, but still, supersymmetry remains unbroken. In section 4, we repeat 
the calculation for the D-dimensional de Sitter spacetime and Kerr-Newman-de Sitter 
spacetimes, finding the same results as in the previous sections. In section 5 we present 
our conclusions along with a brief discussion on some supersymmetry implications for the 
fermion system. 



1 Supersymmetric Quantum Mechanics 

In order to render the article self contained, we briefly review the N = 2 supersymmetric 
quantum mechanics |4 11144 j algebra formalism. For more details we refer the reader to 

inHii]. 

Let H be the self-adjoint Hamiltonian operator describing a quantum system. This quan- 
tum system is supposed to be characterized by the self-adjoint operators {Qi, Qiv}- 
The resulting quantum mechanical system is considered to be supersymmetric, if, 

{Qr.Qj} = H5ij (1) 

with i = 1,2, ...N. The Qi are the corresponding supercharges and the Hamiltonian "iJ" 
is supersymmetric (abbreviated from now on, "susy") Hamiltonian. Therefore, relation 
([1]) describes the N-extended supersymmetry. The anti-commutativity of the supercharges 
results to the following relations, 

2 ^ 

H = 2Ql = 2Ql = ... = 2Q% = -Y.Ql (2) 

i=l 

A supersymmetric quantum system has unbroken supersymmetry, if Eq = 0. However, 
when its ground state does not vanish, that is Eq > 0, susy is broken. 
Making contact with the underlying Hilbert space, susy is unbroken, when 

Q^\i'o) = 0. (3) 

In the above, |^o); denote the ground states. In the special case N = 2, the equation ([2D 
can be written, 

{Qi,Q2} = 0, H = 2Ql = 2Ql = Ql + Ql (4) 
It is more convenient for later purposes, to introduce the operator, 

Q = l={Q^+iQ^) (5) 

and it's adjoint operator, 

Qt = -L(Q,_iQ2) (6) 
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It is easily verifiable that these two satisfy, 



Q2 = Qt2^0, {Q,Q^} = H (7) 

The abstract mathematical theoretical framework of supersymmetric quantum mechanics 
requires an unitary involution defined for the corresponding Hilbert space. We call this 
involution operator, Witten parity and it is an odd, bounded operator [12] which we denote 
W . This operator commutes with the Hamiltonian, 

[W,H]=Q (8) 

and anti-commutes with the supercharges Q and Q\ 

WQ} = {W^,Qt} = o (9) 

Therefore since the Hilbert space corresponding to the Hamiltonian is left invariant by 
the involution W ^ the Hamiltonian H being self adjoint, additionally is an abstract Dirac 
operator. Moreover W serves as an (orthogonal) projection operator on the Hilbert space 
of the quantum system, 

= 1 (10) 

Being a projective operator, W spans the Hilbert space 71 to subspaces classified according 
to their positive and negative Witten parity namely, 

= p^n = ■.w\ij) = ±\ip)} (11) 

Hence, 7i = 7i~^ © T-L~ . Every operator with domain the Hilbert space Ti acquires an 
2N X 2N matrix representation. Using the aforementioned representation we can write: 

with I = icIm, M being the space of N x N matrices. Owing to = and {Q, W} = 0, 
the supercharges can take the following form. 

The matrices A and , are N xN representations of generalized annihilation and creation 
operators, with, 

A-.n^ (14) 

A^ -.n^ ^ %- 



An interesting feature of the representation ()12p . ()13p . is that the Hamiltonian can 
be rendered diagonal, a feature that is rather welcome, due to the fact that the spectral 
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analysis of the Hamilton operator is much simplified. Indeed the spectral problem reduces 
into studying the subspaces [32]. Obviously the Hamiltonian can be written, 

Consequently, the supersymmetric Hamiltonian H, can be written as the direct sum of 
two superpartner Hamiltonians, 

H+ = AA\ H_=A^A (16) 

We shall make use of the operator P^, the eigenstates of which, we call positive and 
negative and we denote |^^). These parity eigenstates satisfy, 

P±|V'±) = ±|V'±) (17) 

and using relation (fl^ . the parity eigenstates can be cast in the form, 

In the above, e H^. Unbroken supersymmetry requires the Hilbert space to contain 
at least one state, \ipo)j with vanishing energy eigenvalue, namely -ff|^o) = 0- Hence, 

Q^m = o (19) 



For a negative parity ground state. 



IV'o>=( ) (20) 



this would obviously imply (by virtue of (jl9p ) that ^I^q ) = 0. Correspondingly, for a 
positive parity ground state, 

m) = ( ' ) (21) 

we would have that ^^|(/>||) = 0. 

A ground state can be either Witten parity positive or Witten parity negative. However 
in a degenerate ground state situation, both cases can occur. 

In the case E ^ 0, the number of positive parity eigenstates equals to the negative parity 
eigenstates. More elaborately, the operator AA^ is unitarily equivalent to the operator 
A^ A, which is equivalent to say that the spectra of the aforementioned operators are equal 
away from zero. This however is not true for the zero modes. In order to describe the net 
number of zero modes, Witten introduced an index, thereafter called the Witten index [31] 
A, 

A = n_ - n+ (22) 
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We denote n± the number of zero modes corresponding to the sub-Hamiltonians H±. 
If the number of zero modes is a finite number, the operator A is Fredholm, which is 
equivalent to saying dimkerA < oo [l2]. Supersymmetry is unbroken if the Witten index 
is a non-zero number G Z. However when the Witten index is zero, we cannot directly 
infer whether supersymmetry is broken or unbroken. The critical issue that determines 
supersymmetry breaking is whether n_|_ = n_ = or if n+ = n„ 7^ |43) . In the first 
case, supersymmetry is broken, while in the second case supersymmetry is unbroken. For 
an interesting application of unbroken supersymmetry when the Witten index is zero see 
references [151 - 147) . 

The Witten index is related to the Fredholm index of the operator A as, 

A = mdA = dimker//_ — dimker//-|_ (23) 

with, 

indyl = dimkerA — dimkervl^ = dimker^^A — dimkerAA^ (24) 

In this paper we shall not make use of Fredholm operators. Therefore we shall need a 
generalization of the Fredholm index (and of the corresponding Witten index). The heat- 
kernel regularized index, both for the operator A, that is indj^ and for the Witten index, 
At, are defined as follows: 

indtA = Tt{-We-'^^^) = tr_(-PFe-*^^^) - tv+{-We-'^^^ ) (25) 
At = lim indt^ 

In the above, t > 0, and additionally the tr-t stands for the trace in the subspaces H^. 
The heat-kernel regularized index is defined for operators that are trace class (in our case 
above, the operator tr(— T^e"*^^"^) must be trace class), that is, they have a finite trace 
norm. This is independent of the orthonormal basis describing the Hilbert space. From the 
Banach space of all trace class operators, we shall be interested in the subspace spanned 
by the A and A^ and their product AA^. 



2 Massless Fermion Quasinormal modes for a Kerr Black 
Hole 

We shall study first the quasinormal modes of a Dirac fermion in a Kerr black hole back- 
ground. Particularly, we focus on the massless fermion case. As we mentioned earlier, 
the Kerr black hole is one of the most important four dimensional black hole solutions, 
since realistic astrophysical black holes are rotating with negligible electric charge and in 
addition, the quasinormal modes stemming from such a background are very important 
for observations of gravitational waves |48fl53j (and more likely gravity waves will come 
from such objects). We shall make use of the Newman- Penrose formalism [53], which is 
frequently used in the literature [25ti32] . Let the metric be (we use the notation of [25]), 

g{ea,eb) = rjab (26) 
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with 



The operators satisfy 



(27) 



(28) 



The Newman-Penrose formalism is based on the tetrads, which is the basis set l,n,m,m. 
The first two are real and the last two complex with the m being the complex conjugate 
of m. These satisfy: 



l-l = n- n = m- m = m- m- 
l-m = n- rh = n- m = n- fh 







(29) 



and I ■ n = 1, m ■ m = 1. The fundamental matrix r]ab in that basis, takes the form: 



ab 



/ 1 \ 

10 

-1 

V -1 / 



(30) 



with the basis ei = /, 62 = n, 63 = m, 64 = m. When used to define directional derivatives, 
these are denoted as: ei = D, 62 = D*, = 6, 64 = 6*. The massless Dirac equation in 
the null tetrad basis reads: 

ij^'V^^ = (31) 
with the covariant derivative being equal to 



v„, = d, 



'm - 7^/. bVcal 



.,cb 



(32) 



The spin connection lo'^^ ^ on the pseudo-Riemannian manifold, satisfies the following equa- 
tion: 

(33) 



For the purposes of our analysis, we write the Dirac spinor ^ in terms of two Weyl spinors, 
^ and ip, namely: 

Substituting Eg . (1341) to Eq . (|31|) . we obtain a new set of equations, in terms of Weyl spinors: 

^^cV^V'" = a^,V^^^ = (35) 

where, represents the Pauli matrices, and and are the covariant derivatives of 
and V'"; respectively. The Pauli matrices take new forms in the vielbein formalism for 
curved spacetime, namely: 



(36) 
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The four dimensional Kerr background spacetime has the following metric: 
, n / 2Mr\ , n /4Mrasm'^9\ , 

ds' = -(l - ^)dt' - ( )dtde (37) 

+ ^dr^ + p'd9^ + (r^ + + '^™'^ ) sin^ 

with, = + cos^ 0, A = + — 2Mr and a = jj the rotation parameter of the 
black hole. For later purposes, we define p = r + ia cos 9 and p* = r — ia cos 9 such that 
= pp* . The basis of the Newman-Penrose formalism, defined for the Kerr geometry can 
be written: 

= ^(A,-p^0,-aAsin2 6'), = ^(A, /o^ 0, -aA sin^ 6^) (38) 

rrij = {—iasm9, 0, —p^, —i{r'^ + a^) sin 6), rhj = m* 

The Ricci rotation coefficients in terms of the spin connection are |54] : 

1 

« = W3ii, p = a;3i4, e = -(a;2ii+a;34i) (39) 

O" = W313, ^ = ^^243, 7=2(^213+^324) 

A = a;244, T = L0^i2i a = -(aJ213+aJ344) 

1 

Using the above, the Dirac equation ()3ip in the Kerr spacetime can be cast in the form: 

{D + e- p)ip° + (5* + vr - a)ip^ = 0, (40) 
{D* + p- 7)^1 + (<5 + /3 - r)^o = 0, 
{D + e* - p*)i}^ + {5 + TT* - Q*)^;^ = 
{D* +p* - + {6* +(3* - T*)tp^ = 



2^ = ^242, 7r = W241, /3 = - (^213+^343 ) 



We can express the basis vectors of the Newman-Penrose formalism as follows: 

A 

1 1 

m = 6 = _ Ci, m = 5* = —pz — £n 
V2p ° V2p* 



l = D = Vo, n = D* = -^vl (41) 



where, 



Vo = dr + i^, Vl = dr-i^ (42) 

Co = dg + Q, Co = dg -Q 
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and K = {r"^ +a?)uj+a'm and Q = auj sm6+-^^. For later use, we introduce the operators: 
^ ^ K r - M . ^ K r - M 

Vn = dr + i-^ + 2n^—, Vi = dr-i— + 2n^^ (43) 

^ „ ^ n ma sin ^+ „ ^ n inasinO 
^"=* + <5+j;^ — 4. = *-'3 + j^ + ^^ 

We seek solution of the form ■ip{t,r,6,(f)) = e'^^^^~^"^'f'^ip{r,9). Upon re-scaling the spinors 
in the following way: 

^{r,e)=p*^p', /2(r,e)=V' (44) 

the Dirac equation can be recast as: 

Voh{r,9) + ^£if2{r,9) = (45) 
AV\f2{r,e)-V2C\fi{r,e) = 

2 2 

Vog2{r,e)-^C\gi{r,9) = 
AV\gi{r,9) + V2Cifi{r,0) = 

2 2 

where Di and Ci, can be deduced from Eq.(|43p. In order to make contact with the 

2 2 

quasinormal modes spectrum, we separate the above functions into radial and angular 
parts, as follows: 

fiir,9)=R_iir)S_ii9), (46) 
f2{r,9)=R,{r)Sii9), 

2 2 

gi{r,e) = R,{r)S_i{9), 
g2ir,9) = R_iir)Si{9) 

By substituting Eq.([lHD, to Eq. (j^S|) . we obtain the following equations: 

£iSi(9) = -XS_i(9) (47) 

2 2 2 

J^\S_i{9) = XSi{9) 



and 



VoR'_^{r) = ^Ri{r) (48) 
vlR.ir) = 'R'_,^ir) 
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with R! 1 (r) = -^R_i{r). Equations (I47p and (I48p correspond to the angular and radial 

~2 VA 2 

part of the Dirac equation respectively. Eliminating Si or S_i from the above equation, 



2 2 



and upon defining u = cos 0, we can find that the angular equation can be written as (with 



, ,A-u^)dSs(6)du] + [iaujuY-2aujsu + s + Aim 5- 5"^ = (49) 

du\ /V 1 — / 

with = + 2maa — {au))"^. The parameter A is real, but our results are indifferent 
to whether A is real or not. This would change slightly our notation, but the results would 
be the same. Equations (j47p and (j48p . will be our starting point of our analysis. We shall 
see that the solutions of these equations are related to a supersymmetric Hilbert space. 
The above equation ([18|) (using the tortoise coordinate defined as = %^)) can be 
transformed to the following Schroedinger like equation: 



d^Z. 



dx^ 



^ + L'^ -V±{x)]Z± = {) (50) 



with Z± = Pi ± P_i and V±{x) = A^-^ it -^l. Equation ([50]) is the quasinormal 



K"^ dx \ K 

modes master equation for a fermion field in Kerr spacetime. We are not interested 
in solving the master equation, there are quite rigorous techniques for doing that [1]- 
I32j . Our main interest is whether the spectrum in terms of uj is discrete or continuous. 
Quasinormal modes are solutions of the above equation, with the wave functions satisfying 
certain boundary conditions at the horizon and at infinity. In addition the quasinormal 
modes corresponding to Kerr black holes form a countable set of discrete frequencies. 
A noticeable feature is the fact that accuracy is necessary when computing quasinormal 
modes because a considerable change in the black hole parameters results in a minor 
change in the quasinormal modes. The boundary conditions are very crucial in order to 
define a trace class operator, which be valuable to us in the following. Based on equations 
()47p and ()48p . we can construct a = 2 supersymmetric quantum algebra acting on the 
fermionic solutions. We define the matrix Dk, 

acting on the vector: 

/ E' , (r\ \ 

(52) 

Replacing the operators Dq from equation Dk is equal to: 

Dk=\ I (53) 
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We can easily obtain Dj^ which is equal to (note that K contains io which is complex): 



■K' \ A 



acting on 




(55) 



Obviously each quasinormal mode satisfies the zero mode equation of Dk- Therefore 
we could say that the quasinormal modes (which can be found from equation ()50p ) are 
in bijective correspondence to the zero modes of the operator Dk- The same argument 
applies for the operator d\^, with the difference that in this case the quasinormal modes 
are the complex conjugate of the previous case. 

Making use of the operators Dk and we can form an = 2 susy QM algebra. 
The supercharges of this algebra, Q and are defined in terms of Dk and Z)]^, 

Dk \ f ^ 



Additionally, the quantum Hamiltonian can be cast in following diagonal form (see previ- 
ous section), 

iZ=f^-^- ) (57) 

The supercharges (|56|) and the Hamiltonian and (|57|) . satisfy the relations d?]) and ([SJ, 
namely {Q, Qt} = h, = 0, Q^^ = 0, {Q, W} = 0, = I and [W, H] = 0. Hence an 
N = 2 susy QM algebra underlies the radial part of fermionic black hole system. 
The equations of the quasinormal modes, Dk\4>~) = and it's conjugate D\^\(f>'^) = have 
complex conjugate solutions. Obviously we can have a bijective correspondence between 
the quasinormal modes, given by equation Dk\4'~) = 0, and their complex conjugate 
counterparts, given by D^j^\(f)'^) = 0. It is obvious that this bijective correspondence holds 
between the zero modes of the matrices Dk and d\^. Therefore, 

kerL>i^ = kevD]^ (58) 

which in turn implies kerD^^D^^ = keiDxD^x- ^ consequence, the following relation 
holds for the operators e~^^K^'^ and e^*-^^^^ 



= tr+e-*^^^!^ (59) 



tr_e" 

Recall that tr-t stands for the trace in the subspaces H^. As a consequence of relation ()59p . 
the regularized index of the operator Dk is equal to zero. Consequently the regularized 
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Witten index is also zero. Hence, the fermionic system possesses an unbroken N = 2 
supersymmetric quantum mechanics algebra. 

Using the notation of relation ()20p the equation = has a direct representative 

equation for the supercharge, namely, Q\iPq) = 0. In view of section's 1 terminology, this 
implies that the zero mode eigenstate IV'q ) is a negative Witten parity eigenstate, and is 
equal to: 






R 



\ 



(60) 



In the same vain, the positive Witten parity vacuum eigenstate is. 



/ (r' ,1 



i?,i(r 



V 



\ 



(61) 



The angular case can be treated accordingly following the same lines as in the radial part 
of the Dirac equation. Adopting the notation of the previous paragraphs, the algebra can 
be built on the matrices Djig and D^^g- "^^^ matrix D^g is defined to be: 



Ci 



(62) 



acting on 



8.(9) 



(63) 



which is a direct consequence of equation ()T7|) . It's conjugate equals to: 



ok 



£1 A 
-A CV 



(64) 



Recall that the operators £i and C\ contain ui which is a complex number. We must 

2 2 

note that the situation at hand is much more evolved in comparison to the angular case. 
It is obvious that since quasinormal modes exist, the operator Dxg certainly has a set 
of discrete zero modes, that belong to a countable set of complex numbers. However we 
cannot argue that the same holds for the operator Z)]^^ . In the radial case we could solve 
the zero mode problem of the two corresponding operators simultaneously, since the zero 
modes of the operators where complex conjugate to each other (a proof for existence of 
solutions in a much more general setup see [59]). But in the angular case we cannot apply 
the same argument. Therefore we can argue that in general the number of zero modes 



12 



of the two matrices are not equal. Hence we could naively argue that super symmetry is 
unbroken in this case, but for different reasons in comparison to the radial case. This 
naive argument is based on the fact that the operators have not the same number of zero 
modes, and therefore supersymmetry is unbroken. This however would be true only in 
the case the operators where Fredholm, which are not (since dimkerD^g — )• oo). Moreover, 
we cannot be sure whether the operator dJ^^ is trace-class. Therefore, we conclude that 
only the radial part of a Dirac fermionic system in the Kerr black hole background can be 
associated to a = 2 susy QM algebra. 

Before closing we must note that up to date, the important theoretical issue that addresses 
the nature of neutrino, that is whether it is Dirac or Majorana, has not be answered 
successfully yet. Hence any information on the effect of massless (if the neutrino can 
be considered massless) fermions in nature is invaluable. We studied a massless Dirac 
fermion in the most realistic curved four dimensional gravitational background and found 
an underlying supersymmetry, a result that has some value in combination that is very 
closely related to the quasinormal modes of the fermion in the same background. It would 
be certainly be interesting to study if super symmetric structures exist when massive Dirac 
fermions and also when Majorana fermions are considered. Of course this would require 
a complete study of the quasinormal modes of massive Dirac fermions and of Majorana 
fermions in Kerr backgrounds. We hope to address such issues in the future, but these 
issues exceed the purposes of this paper. 



3 Fermion Quasinormal modes for Kerr-Newman, Reisnerr- 
Nordstrom and Schwarzschild Black Holes 

3.1 The Kerr-Newman Black Hole 

In this section we further explore whether supersymmetric structures underlie any other 
fermion systems in curved gravitational backgrounds. We shall study first the Kerr- 
Newman black hole, which is the only asymptotically flat solution of the Einstein equa- 
tions, with electrised vacuum. We consider a massive fermion of Dirac type in such a black 
hole background. Following the steps of the previous section and adopting the notation 
of reference [26], the null tetrad corresponding to Kerr-Newman black hole are: 

I'' = ^{r^ + a^,A,0,a), n'^ = -^(r^ + a^, -A, 0, a) (65) 
A 2pp 

1 i 
= ^ (iasinO, 0, 1, 

with, p = r — ia cos 9, IS. = r"^ + — 2Mr -\- Q^. Using as previously the Newman-Penrose 
formalism, the general Dirac equation for massive fermion in a Kerr-Newmann background 
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reads, 



{D + e- 


- p)P° + {5 + TT- 


a)P^ -- 




(A + /U 


-j)P^ + {5 + /3- 






(D + e- 


- p)Q^' + {5 + ^- 




= —impP^ 

V2 


(A + /X 


-7)Q'' + (^ + /3- 


- f)Q'' 


= -^impP^ 
V2 



(66) 



We write the wave functions in a form that the radial and the angular part are separated. 
Then the wave functions can be written as follows, 

P° = -R_ 1 {r)S_ 1 {9) exp-^("'*-'^<^) (67) 

p 2 2 

= ip, i(r)5, i(^)exp-^(^*-™<^) 

p ' 2 '2 
P ^ 2 2 

Q°' = --R_i{r)S.i{e) exp-^("*-"^<^) 

P 2 ^2 

Hence, using the equations ([66]) . Dirac can be recast as, 

CiSi = —{\ — amp cos9)S_i (68) 

2 2 2 

C\ S_i = (A + amp COS 6)Si 

2 2 2 

in reference to the angular part. The radial part can be written as: 

VAVoR_i = (A + impr)R!. (69) 

2 +2 

VapJp'^i = (A - impr)V2R^i 
with = a/AP^i The operators that appear in relations (f68|) and ([69|) are equal to: 

2 2 

_ ^ K n dA _x ^ if n dA 
^'^ = ^^ + ^A + Adr' ^^ = ^^-^A + Ad7 

Cn = do — aui sin -| + n cot 0, C}, = de + au; sin 6* + n cot 

sm 8 sm 

with if = (r'^ + 0^)0; + am. Note that A is the same as in the Kerr case, since the angular 
equation can be reduced to Eq. (fl9|) . corresponding to the Kerr case [26]. In the following 
paragraphs we shall study both the massless and massive fermion case. As we shall see, 
the mass can introduce some complications to our initial arguments but it remains robust 
to these complications. 
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3.1.1 Massless Fermion case 



Consider the radial part of a massless Dirac fermion first. As in the previous section, from 
the equations of motion (j69p , we can construct the matrix Dpt on which the supersymmetric 
quantum algebra can be built on. This matrix is defined to be: 



Vadq a 



Dr={ ' , ' ^^^t 1 (71) 



acting on the vector: 



(72) 



/ R_.{r) \ 

Using the explicit form of the operators defined in equation ([70|). the operator Dr equals 
to: 

We can easily obtain it's adjoint, L*^, which is equal to: 



acting on 



(75) 




The number of the zero modes of the operator are bijectively related to the number of 
the zero modes of D^. This is because, the set of the zero modes of are in one-to-one 
correspondence to the quasinormal modes, corresponding to the equation ([69]) . In the 
same vain, the zero modes of are in one-to-one correspondence to the quasinormal 
modes, corresponding to complex conjugate of the equation (|69|) . It is necessary to note in 
order to obtain consistent solutions for the complex conjugate of equation (|69p . the wave 
functions must obey the complex conjugate boundary conditions of the wave functions 
corresponding to equation ([69]) . The existence of a solution for this case is obvious but 
can be further justified by a theorem on second order differential equations [59]. We are 
not interested in the specific form of the solutions, but only in the fact that the zero 
modes of the operators and Dji have a bijective correspondence. The situation at 
hand is very similar to the massless Kerr fermion case of the previous section. Having 
found a correspondence between the zero modes, it is easy to verify that supersymmetry 
is unbroken, with the heat-kernel regularized Witten index being equal to zero. Let us see 
this in detail. The susy QM algebra can be built on the supercharges, 
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and also the corresponding Hamiltonian, 

These satisfy {Q,Q^ = H, = 0, Qt^ = q, {Q,W} = and [W,H] = 0. Hence an 
N = 2 susy QM algebra underlies the radial part of fermionic Kerr-Newman black hole 
system. In addition, as in the previous section, as a result of the bijective correspondence 
of the two matrices zero modes we have, 

keiDn = kerL>^ (78) 

which in turn implies, kerDfiD^j = kerD/jD^j. As a consequence, the following relation 
holds for the operators e~^^R^^ and e~^^^^R 

tr_e-*^«^« = tr+e-^^^^R (79) 

The above result holds only in the case the operator ti^—We^^^R) is trace class, which is 
true, since the operators Dj^, and consequently the operators DjiD\^, D^^Db. are trace 
class [[42] . Hence the regularized Witten index is zero and this in conjunction with equation 
()78p implies that the radial part of the fermionic system has unbroken supersymmetry. 
We omit the study of the angular part of the fermionic system since as in the Kerr case, 
supersymmetry is broken, for the same reasons. 



3.1.2 A brief Discussion 



We mentioned in the above that in order relation ()79p holds, the operator tr(— VFe^^^fl) 
must be trace class. We argued that this is true since and Z)Jj are trace class. 

However this is not prerequisite for our case, and the only requirement for tr(— H^e^^^fl) 
to be trace class is that Dr L'Jj — Dr is trace class. Indeed, by virtue of a theorem [12] , 
if T and S are two self-adjoint operators, and T — 5 is trace class, then /(T) — f{S) is 
also trace class. The function / is a map f : TZ ^ TZ, satisfying: 

POO 

/ /(P)(l + |p|)dp <oo (80) 

J oo 

with / the fourier transform of /. It is clear that in our case, the function / is the trace of 
the exponential of the operators Dji Z)Jj and Dr. Hence the theorem applies perfectly 
in the present situation. 

3.1.3 Massive Fermion case 

We extend the study in the case the Dirac fermion has mass mp. From the equations of 
motion ()69p we can construct the matrix, 

y/A(dr + i^) X + impr 

X-irrJ M^r-il' I '''' 
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which acts on the vector: 

The adjoint of the matrix Dmr is equal to: 



(82) 



\/Aydr — i^j \ + impr 
X — impr \^(dr + i^ 



Dln= I I (83) 



Unhke the massless case, it is not obvious that the zero modes of exist at all. 

One must solve the equation D^jj^ip = 0, subject to the complex conjugate boundary 
conditions corresponding to equation ()69p . The zero modes of the matrix Dmr exist and 
the countable set that these belong to is in one-to-one correspondence to the quasinormal 
modes for a massive fermion in the Kerr-Newman background. The latter were studied 
in detail in [26]. We are thus confronted with a non trivial problem. In any case, there 
are two alternative situations that can occur, in reference to the zero modes of d\,^^. 
Either dimkerZ)]^^^ = or dimkerZ)]^^^ 7^ and we cannot be sure which case is true, 
unless we solve explicitly the equation D'^^^^ijj = 0. Nonetheless we can answer whether 
supersymmetry is broken or not if we make use of a theorem. First note that we can write 
Dmr = Dn + C, with C the odd symmetric matrix |42j . 

C=( .0 ^"^/^"j (84) 

and also, the Dn is the massless Kerr-Newman case matrix defined in equation ()73p . Now, 
if the operator iiW e"^^^^^'' is trace class, the following theorem holds, 

indt{D + C) = \ndtD (85) 

with C a symmetric odd operator and D any trace class operator. That is, the regularized 
index of the operator D + C \s equal to the index of the operator D. In our case, this 
reads, vadt{Dpi + C) = indfD/j. Recall that we found earlier in this section that the Witten 
index and the heat kernel regularized index of D^j in the Kerr-Newman background are 
both zero, hence indi-D/j = 0. By virtue of theorem (j85|) . indt(Z)ij + C") = 0. Hence, in 
conjunction with the fact that kerZ)/? 7^ 0, we may argue that supersymmetry is unbroken 
even for the massive fermion case. 

Before closing we must mention that if we could be able to prove that dimkerDj^^^^ = 0, 
then supersymmetry would be unbroken. However we gave a different proof that an 
unbroken N = 2 supersymmetric quantum mechanics algebra underlies the Kerr-Newman 
fermionic system. 



3.2 The Reisnerr-Nordstrom and Schwarzschild Black Hole 

Following the lines of argument of the previous sections, it can be easily proven that an 
unbroken N = 2 supersymmetric quantum mechanics algebra underlies the radial part of 
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massless and massive Dirac fermionic systems in Reisnerr-Nordstrom and Schwarzschild 
black hole backgrounds. For a detailed study of the corresponding fermionic quasinormal 
modes can be found in references [27] and p8], respectively. We shall not get into details, 
since the results are identical to those of the Kerr-Newman case. Indeed, the radial part of 
the Dirac equation in the Reisnerr-Nordstrom background can be reduced to the same form 
of equations, like Equation (j69p . In the Reisnerr-Nordstrom black hole, A = r"^ — 2Mr+Q'^ 
and = + with j positive. The charges of the super symmetric quantum mechanical 
algebra are identical to those of the Kerr-Newmann case, with the appropriate replacement 
of A and A. Definitely, the boundary conditions for the functions to yield the corresponding 
quasinormal modes must be taken into account, but the reasoning is the same. In the case 
of the Schwarzschild black hole the same arguments hold, but we should replace everywhere 
A = — 2Mr and A^ = (j + |)^. Before closing this section we should mention that 
someone could argue that the supersymmetric quantum mechanics algebras corresponding 
to the Reisnerr-Nordstrom, Schwarzschild and Kerr-Newman problems, could be related 
by virtue of the fact that the Reisnerr-Nordstrom metric is the Kerr-Newmann black one 
when a = 0. In addition, the Schwarzschild black hole is the is the Kerr-Newmann black 
hole when a = and (5 = 0. This would be a very interesting task, but we defer such an 
issue to a future work. 



4 The D-dimensional de Sitter spacetime and Kerr-Newman- 
de Sitter Black Hole 



In this section we shall study whether the supersymmetric structure we found in the 
previous sections, also underlies de Sitter related spacetimes. We start first with a D- 
dimensional de Sitter spacetime, in which we consider a massive Dirac fermion. We adopt 
the notation of reference [29]. See also references [30H32j . In general, the metric of an 
D-dimensional spherically symmetric spacetime, is of the form, 



ds^ = W'^dt^ 



1^ . 



D-2 



(86) 



with dn^_2> the metric of the D — 2 dimensional unit sphere. In this spacetime, the Dirac 
equation of a fermion with mass ruf can be cast as: 



rufW 

.TTd_,t.V 



^4 + ^7 

— TUfW 



( fi%) 



(2), 



UJ 



(87) 



with k = ±(— 2 1- n). Now in the D-dimensional de Sitter spacetime background, the 



quantities U, V and W are equal to: 

V _ 1 
y Lsinh(^) 



U= 1 



w 



Lcosh{M-) 
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In the above equation, L is related to the cosmological constant. Using the coordinate 
z = tanh(^), with r = zL, the Dirac equation (j87p can be written: 



lOJ 



dz ^ 
d 



lUJ 



Az ^ 



Ri{y) 
R2{y) 



k 



iz 



iM 



[k 

- +iM 

z 



i?2(y) = 
^i(y) = 



(89) 



with uj = u) L and M = nif L and additionally 

Riiy) = -ifiiiy) + fi%)^ My) 

In the case M = 0, equation ([89]) becomes: 



^ £ky) - f:.'.L(y) 



(2), 



vi-^ 



dz 71^^ 



k 



d 



dz ^ 



Ri{y) - -R2{y) = 
R2{y)--Ri{y) = 



(90) 



(91) 



The above equation gives the quasinormal mode spectrum of the fermionic system in this 
de-Sitter background. It is obvious that also in this case, as in the previous cases, the 
quasinormal spectrum is in bijective correspondence to the zero modes of the operator: 



Dds 

acting on the vector: 
and also D^^^g- 



y2_d_ iu 



VT 



z 

Z2 d 



VT- 



Ri{y) 
R2{y) 



(92) 



(93) 



D 



^/^'lz + 



dS 



VT- 



acting on 



R2{y) 
Ri{y) 



(94) 



(95) 



As in the Kerr case, the operators can be used to form an unbroken N = 2 supersymmetric 
quantum mechanics algebra. Without getting into details, the supercharges of this algebra, 
Q and are. 



Q 



Dds 








(96) 
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Moreover, the quantum Hamiltonian is. 



H = 




) 



(97) 



Since 



kevD 



dS 



(98) 



which in turn imphes, kerDdsD'^^ = keiDdsD'^g. As a consequence, the fohowing relation 
holds for the operators e~^^ds^'is and e~^^'^^^ds 



Hence supersymmetry is unbroken, just in the Kerr massless fermion case. The difference 
between the two problems however is that in the de Sitter case, the whole system possesses 
this supersymmetry and in the Kerr case, only the radial part has this symmetry. The 
massive case can be also treated identically to the massive case of the Kerr-Newman black 
hole, yielding the result that the system possesses unbroken N = 2 susy QM. Before closing 
this section we mention that the system of a Dirac fermionic field in the Kerr-Newman-de 
Sitter black hole background has also an = 2 susy QM, but we omit such a study since 
it is identical to the other cases we studied, with different supercharges and Hamiltonian 
of course. For a study of the fermionic quasinormal modes in a Kerr-Newman-de Sitter 
background, see reference [30] . 

5 Discussion and Conclusions 

In this paper we studied the relation of the quasinormal modes to supersymmetry. In par- 
ticular, we showed that an = 2 supersymmetric quantum mechanics algebra underlies 
Dirac fermionic systems when these are considered in various black holes backgrounds. We 
examined thoroughly the massless Dirac fermion Kerr black hole, the massless and massive 
fermion case in Kerr-Newman, Reisnerr-Nordstrom, Schwarzschild, D-dimensional de Sit- 
ter, and Kerr-Newman-de Sitter spacetimes. In all the cases we found a hidden unbroken 
supersymmetry in the radial part of the Dirac fermion system, and the non-breaking of 
supersymmetry was very closely connected to the operators being trace-class and to the 
existence of quasinormal modes. Actually the zero modes of the fermionic system in each 
case, was in bijective correspondence to the quasinormal modes of each system. 
Let us see some implications of supersymmetry for the fermionic systems. Firstly, from a 
mathematical point of view, some quantities related to the corresponding spectral problems 
have a very simple form, for example the Krein's spectral function ^(A) |42j . which is zero 
since the regularized Witten index is zero. Additionally, the Hamiltonian is invariant 
under the global U{1) transformation of the supercharges [43j. 



Therefore, the resulting quantum system possesses a continuous R-symmetry, the breaking 
of which into a discrete one can be a very interesting situation. An intriguing question is. 



(99) 



Q' 



(100) 
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how to break this continuous symmetry exphcitly, because an exphcit breaking could make 
us assume that a supercharge acquires a constant vacuum expectation value (assuming it's 
vacuum expectation value is it's trace over all vacuum eigenstates). Such a thing does not 
spoil the trace-class properties of the operators we studied and it's physically appealing 
since discrete symmetries are inherent to quantum system with fermionic condensates. 
In view of gauge/gravity dualities it would be interesting to search for such supersym- 
metries in quantum systems with AdS backgrounds. In addition, such a study would be 
complete if both Majorana and Dirac fermion zero modes are studied in such backgrounds, 
in reference to supersymmetric structures. Moreover this would require the Majorana 
fermion quasinormal mode spectrum, in black hole backgrounds, and since the nature of 
the neutrino is not known (that is whether it is Majorana or Dirac), up to date, it would 
be interesting by it self to study this issue. In addition, classical black holes solutions may 
be used in the quantum black hole theory. A perturbative treatment of these classical 
solutions is a first approach to a quantum black hole [60]. It is well known [60] that non- 
Abelian uplifted magnetic dilatonic black holes are connected to a = 4 supersymmetric 
quantum mechanics algebra. A question that naturally arises is whether this A^ = 4 al- 
gebra can in any way be connected to the N = 2 algebra we studied in this article. It 
is obvious that the theoretical frameworks are different in the two cases, but nevertheless 
the underlying supersymmetric quantum mechanics structure motivates us to check also 
whether there is a connection between the topology of the black holes and the N = 2 susy 
QM algebra we found. This connection to the topology is also advocated by the fact that 
the Witten index of the algebra vanishes, the last being a topological property character- 
izing the various spin structures of the spacetime. We hope to address these issues in a 
future work. 
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